The Evolution of Inverse Power Law Quintessence at Low Redshift 
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Quintessence models based on a scalar field, <f>, with an inverse power law potential display simple 
tracking behavior at early times, when the quintessence energy density, p^, is sub-dominant. At late 
times, when p^ becomes comparable to the matter density, p m , the evolution of <f> diverges from its 
scaling behavior. We calculate the first order departure of (j> from its tracker solution at low redshift. 
Our results for the evolution of <f>, p^, f2^, and w are surprisingly accurate even down to z = 0. We 
find that w and Q$ are related linearly to first order, and derive a semi-analytic expression for w(z) 
which is accurate to within a few percent. Our analytic techniques are potentially applicable to any 
quintessence model in which the quintessence component comes to dominate at late times. 



PACS numbers: 98.80.Cq 



I. INTRODUCTION 

A host of observations suggest that the universe is cur- 
rently dominated by some form of energy with negative 
pressure El 0> Hi • Recent observations of Type la su- 
pernovae jj| , which appear fainter than they would if the 
expansion were decelerating due to matter alone, provide 
the most direct evidence for a dominant negative pressure 
component (and the accelerating expansion that accom- 
panies it). In addition, cosmic microwave background 
observations from WMAP, for example, suggest that the 
universe is nearly flat (Oo re 1) but contains a sub-critical 
density of matter Q m o = 0.27 ± 0.04 0- Taken together 
these findings imply the presence of a large fraction of 
unclustered "dark energy" . 

The most straightforward source of dark energy is a 
cosmological constant, but the level of fine-tuning nec- 
essary to make it a viable candidate is disconcerting. 
An alternative, dubbed quintessence, is a model of time- 
varying dark energy based on the behavior of a dynami- 
cal, classical scalar field <j>. In the simplest quintessence 
models, <j> is taken to be a minimally-coupled field with 
potential V(<j>), and pressure, energy density, and equa- 
tion of state 
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P4> = -^ + V{4>), 



and 



w 



(1) 



(2) 



(3) 



The behavior of the quintessence field will depend, of 
course, on the particular choice of V(</>). One widely- 
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investigated case is the inverse power law potential 0] 



V{4>) = K<t>~ 



(4) 



-,4+a 



where k is a constant with units of rrr 

These potentials are particularly simple and have sev- 
eral desirable properties. They exhibit tracking behav- 
ior, i.e., a wide range of initial conditions converge to a 
common solution at late times 0, 0- I n addition, in- 
verse power law quintessence remains sufficiently sub- 
dominant energetically so that it does not disastrously 
interfere with known epochs of standard big bang cos- 
mology. Finally, because a fundamental understanding 
of the nature of dark energy must eventually come from 
particle physics, it is worth noting that inverse power 
law potentials can arise naturally from extensions to the 
standard model (e.g., SUSY QCD Q). 

The inverse power law tracker solutions are character- 
ized by an equation of state parameter, w, which is ap- 
proximately constant when the universe is dominated by 
a background fluid. During the matter-dominated era, 
for example, when the contribution of the quintessence 
energy density to the expansion is neglected, w is given 
bylEiQ: 



a 



and 6 evolves as 



corresponding to a density which evolves as 



(5) 



(6) 



(7) 



These expressions provide an excellent approximation to 
the behavior of the quintessence field as long as p$ / p m <C 
1. 

At late times, however, when the scalar field itself con- 
tributes significantly to the expansion of the universe, 
the value of w begins to diverge from its tracker value, as 
do <fi(t) and p^(t). In this paper, we examine analytically 
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the behavior of the scalar field during the transition from 
matter to scalar-field domination. It is this late-time evo- 
lution which is relevant to supernova la observations and 
other tests of dark energy (see, e.g., Ref. and ref- 
erences therein). A good analytic approximation to the 
behavior of £1$ and w at late times would be useful in 
such calculations. Furthermore, our treatment provides 
insight into the behavior of the scalar field once it begins 
to dominate. 

A number of observations provide stringent upper lim- 
its on w. For example, the WMAP results, combined 
with supernovae observations, HST data, 2dFGRS obser- 
vations of large scale structure, and Lyman alpha forest 
data give w < —0.78 at 2 a @|, assuming constant w, 
which implies a low value for a. Nonetheless, we have 
attempted to keep our approach as general as possible, 
treating all values of a in the discussion below. 

In the following section, we rederive the inverse power 
law model solutions for 0, p^, and w during the matter- 
dominated era. In Section UTTl we use the first-order per- 
turbations to these background solutions to describe the 
deviation of the scalar field away from tracking behavior 
as the universe enters the scalar-field dominated era. In 
Section HV1 we compare our analytical results to those of 
the numerically-integrated evolution equations. We find 
that the first-order corrections to the background solu- 
tions for 0, p^, f^, and w characterize their behavior 
surprisingly well even to the present. In Section we 
discuss some consequences of these results. 



II. INVERSE POWER LAW QUINTESSENCE 
EVOLUTION IN THE MATTER-DOMINATED 
ERA 



The equation of motion for </> is 



(8) 



where H = a/a is the Hubble parameter and a = a(t) 
is the scale factor which describes the expansion of the 
universe as a function of time. 

We confine our analysis to z < 30, so that the contri- 
bution of radiation to the expansion can be neglected. 
Assuming a flat universe containing only matter and 
quintessence, the Friedmann equation is 



m 2 = Pv 



(9) 



in units where 8nG = 1, which we will use throughout 
the paper. 

When is neglected in equations © and ©, it is 
easy to integrate these equations directly to derive </f>(i), 
as was done in Refs. [EIEI3- However, in deriving a per- 
turbative expansion when p^ / p m is small but nonzero, it 
turns out to be easier to use a as the independent vari- 
able, rather than t. Further, we would like to express 
observable quantities such as w and at late times as 



functions of redshift, rather than time, which is straight- 
forward when the independent variable is taken to be a. 

Appealing to equation (J2J) and equation Q, \<j) 2 can 
be rewritten as 



1 „ 



■{Pm + V), 



2 T 1-x' 
where V implicitly means V(<f>) and 



2( P rn + P4,) 3H 2 



l\ a Ta) 



(10) 



(11) 



Substituting equation (^UJ into Eqs. (JTJ, J2J), ©, and 

©, we can rewrite p^, w, and the Friedmann equation 
as follows: 



P<t> 



xp n 



V 



1-x 



and 



Xp m ~ V(l - 2x) 
Xp m + V 

2 _ Pin + V 



1-x 



(12) 



(13) 



(14) 



Using the above expressions, equation JSJ) becomes 

a<t> + — {5 - 3x) + ^Tv\— + m)=°> (15) 



where the prime denotes d/da. During the matter- 
dominated epoch, when x <C 1, equation (|15fl reduces 
to 



%JJ , , W (0) 3 dV 



>(0) + 



2 p m d<p 



(16) 



where the zero subscript in parentheses denotes the 
zeroth-order solution, neglecting the contribution of 
quintessence to the expansion rate. (We use the subscript 
"(0)" throughout to denote such zeroth-order terms, 
while the subscript "0" without parentheses is reserved 
for quantities evaluated at the present day). 

For quintessence with an inverse power law potential 
(equation 0J, the solution to equation l(T^|) is 



(o) =C(a)a 3 /( 2 +«). 



(17) 



In equation (|17|l . the constant C(a) depends on the 
value of a and is related to n and the present matter 
density, p m o, m the following way 



3(4 + a) , 2+a 

K = p m0 - r^C (a) 

2a(2 + a) z 

The corresponding zeroth-order solution for p^ is 



C{af 



Pmo) — 3pmo /r> ; — r 
y ' a(2 + a) 



-3a/(2+a) 



(18) 



(19) 
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The zeroth-order behavior of the quintessence density pa- 
rameter is then given by Q^ ) = P<t>(p) / (p<t>(o) + Pm)- Fi- 
nally, the equation of state is just a constant in the track- 
ing regime: 



(20) 



These expressions for 0( O )(a), w^(a), and p^(o)i a ) are i 
of course, equivalent to the previously-derived time- 
dependent quantities given in Section I. 



III. FIRST-ORDER PERTURBATION TO 
INVERSE POWER LAW QUINTESSENCE 
EVOLUTION 

We define the first-order perturbation to the zeroth- 
order background field 0( O ) (equation I17fl to be 0(i), 
where 

0(i)/0(o) ~ O(p^ )/p m ). (21) 
Then the perturbed field, 0, is, to first order, 

= 0(o)+0(i): ( 22 ) 

where we use a tilde throughout to denote quantities ex- 
panded to first order in p^fm/pm- Substituting equa- 
tion (|22|l into equation l|15l) and keeping all terms of or- 
der p^ )/Pm (note, in particular, that x ~ O(p^ )/p m )), 
we find that 0^) obeys 



a 0(i) + o( 5 0(i) ~ 3a; (o)0(o)) + 



2 y T ^> '"' r w 2/3. 
3/00(0) 3 d 2 F 



(/5m) 2 d P" 



i) - 0, (23) 



where x (0) = (a0' (Q) ) 2 /6. 

For our particular choice of potential, equation 123|) 
becomes 



«V(i) + ^( 5 0(i) -3^(0) 0(o)) + 

0. 
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3a(l + a)/ta 0m 



/W0 ( o)) (2+tt) 
The solution to equation l|24|l is 



3(6 + a)0 3 



(0) 



(24) 



(25) 



rw a(2 + a)(28 + 8a + a 2 )' 
The corresponding corrections to pMo) and uv ) are 

0(1) 

P0(1) = ^(0)^(0) + I(l)/>m - "7 V, (26) 

0(0) 



and 



^(l) = (1-^(0)) + a l— )> ( 27 ) 

V P0(O) 0(0) y 



where i(x) = 2£(o)0'( 1 )/0'( O ) is defined via a; = iE(q) + 

^(D =« 2 (0fo) +20' ( o) 0(i) )/6. 

In terms of 0( O ) an( i 0(1) j ^ ne first-order perturbations 
to /O0(o) an d ^(o) are 

_ a(4 + q) /0(i)\ . . 

p^-—QT^\j^) Pm ' (28) 



and 



a(A + a) / 0(i) \ 

6 + Q \0(o)/ 



= ^T^-l^J^o)- ( 29 ) 



Because 0m is negative, equation and equation 
predict that pMi) is positive and wn\ is negative, as 
expected. Note also that w^/w^ — p^m/p^fo) — 
P^(i)/2p0(o)- 

Combining the zeroth-order and first-order solutions 
for 0, ptf,, and u>, we get 



<P(p) 

0(0) 



3(6 + a) 3 
a(2 + a)(28 + 8a + a 2 ) <P (°'' 



(30) 



a(2 + a)(28 



: 28 + 8a + a 2 V y /' 



P<j> — P<f>(o) + P«£(l)' 



_ a(4 + a) / 0(i) 
6 + a V0 (o) 



(31) 



n A -L 3(4 + a) r t/^2 n 6/(2+a) N > 

" P ^ I 1 + (2 + a )(28 + 8a + a 2 ) C(a) ° J' 

and 

w = W( ) + 

-->(' 

-- ™( 1 + ptJ^wj c <"> 2 " 6,(2+ " 1 )- 

We also define = p^/{p,f, + p m )- 



IV. A COMPARISON OF ANALYTIC AND 
NUMERICAL RESULTS 

To test the accuracy of our analytical predictions, we 
numerically integrate equation (jSJ) for the cases a = 1, 
a = 2, and a = 4. In order to compare analytic and 
numerical results for a specific case, we adjust k to give 
Q m o = 0.3 and f^o = 0.7. In figures 1 — 4 below, the 
zeroth-order (tracker) and first-order analytic solutions 
for 0, p^, Q<k } and w are compared to their numerical 
counterparts. 

In Fig. 1, we present the ratios of the tracker and first- 
order quintessence fields to the true (numerical) field: 
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</>(o) /<?!>, and <j>/</), respectively. As expected, <fi agrees ex- 
ceptionally well with the true evolution at early times, 
when the quintessence density is sub-dominant; the error 
for all three cases is less than 2% for z > 1. The analytic 
expression begins to break down at later times, becoming 
significantly less accurate at z — 0, but it is still much 
more representative of the true behavior than the tracker 
solution 4>(q). 
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FIG. 1: The ratios 4>(o)/4> (upper three curves) and 
(lower three curves) for a = 1 (black, solid), a — 2 (red 
dotted), and a = 4 (blue, dashed). 



In Fig. 2, we present p^o)/ p<j>, and py,/p0. Here the 
power of the first-order approach is evident. The first- 
order solutions in these cases are accurate to within 6% 
all of the way down to z = 0. In contrast, the tracker 
expression for p$ grossly underestimates the quintessence 
density at low redshift. 

In Fig. 3, we give the ratios fi^ftyjA^ an d ^<>/^0- 
The agreement between Q,^ and the true is quite re- 
markable. Down to z = 1, the error in the first-order 
expression is less than 1% for all three cases we have 
examined, and even at z = 0, the error is less than 2%. 

Finally, we consider the equation of state parameter, 
w. In Fig. 4, we show the ratio of the tracker value 
(which is, of course, constant) to the true w and w/w. 
Agreement between w and w is quite good at early times, 
with an error of less than 2% for z > 1, increasing to 4% 
at z = 0.5 and 10% at z = 0. 



V. DISCUSSION 

Our first-order corrections to the pure matter- 
dominated solutions for <t>, p&, CIa, and w characterize the 
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FIG. 2: The ratios p^,(o)/p<t> (lower three curves) and p^fp^ 
(upper three curves) for a = 1 (black, solid), a = 2 (red, 
dotted), and a = 4 (blue, dashed). 



evolution of these quantities extremely well at 1 + z > 2. 
Although, unsurprisingly, the first-order solutions begin 
to break down for 1 + z < 2, when the quintessence com- 
ponent dominates the expansion, they are suprisingly ac- 
curate even to the present. 

We can draw some simple quantitative conclusions 
from these results. It is instructive to use equations H19|) 
and l|32f> to express w as a function of f2^, since these 
are the two directly-observable quantities of interest. We 
find, to first order in p<j,(a)/pm' 

2a(4 + a) 



w = — ?rr"ii- (33) 

2 + a (2 + a)(28 + 8a + a 2 ) V ' 

The first term in this equation is just a)™, while the 
second term shows how w diverges from the tracker value 
at early times, when the quintessence field is just starting 
to contribute to the expansion of the universe. We see 
that the deviation of w from W( ) depends linearly on 
Further, the constant of proportionality is a slowly- 
varying function of a, ranging from 0.09 for a = 1 to 0.14 
for a = 4. Hence, the rate at which w deviates from w/q\ 
at early times is almost independent of a. Equation 133|) 
is an excellent fit for small O^, but it works suprisingly 
well even up to the present. For a = 1 — 4, the error in 
the prediction for w is less than 10%, with smaller errors 
for smaller (and more physically relevant) values of a. 

A number of approximations have been proposed to 
simulate the evolution of w at late times. For arbi- 
trary quintessence models, a common app roximation is a 
Taylor expansion in z 



lcis, a common appr 

nam mm nt: 



Win 



W + WlZ + 



(34) 
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FIG. 3: The ratios n^oj/fi^ (lower three curves) and fi^/n^ 
(upper three curves) for a — 1 (black, solid), a = 2 (red, 
dotted), and a — 4 (blue, dashed). 
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FIG. 4: The ratios W(o)/w (lower three curves) and w/w 
(upper three curves) for a — 1 (black, solid), a — 2 (red, 
dotted), and a — 4 (blue, dashed). 



while more complex parametrizations have been exam- 
ined by Bassett et al. [jja] and by Corasaniti and Copeland 

M 

For the specific case of inverse power-law models, Efs- 
tathiou |3 showed that a good empirical fit to the time- 



WEfst = w +wi ln(l + z). 



(35) 



Our results suggest an alternative parametrization. We 
can write equation i|32|l in the form 



wet = «; (0 ) + - «; (0) )(1 + z)- 6 /( 2+Q ), 



(36) 



Since u>( ) is given by equation 1)21) [). equation (|36[l has 
only one undetermined parameter, the present value of 
the equation of state, wq, which can be determined nu- 
merically and used to derive a fit to w. As shown in Fig. 
5, the error introduced by this approximation is less than 
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The ratio wat/w, where wa t is given by equation 
v — i fwiar-Tr ar.KA"\ n — 2 (red, dotted), and a = 4 



FIG. 5 

I36fl , for a — 1 (black, solid), a 
(blue, dashed) 



3% for a in the range 1 — 4. 

To show how our fit (eauation l36(l compares to the em- 
pirical fit of equation (|35|) , we have plotted their behavior 
for the a = 1 case in Fig. 6. A x 2 minimization routine 
was used to determine the best-fit value of W\ over the 
plotted redshift range for u>Efst- (In both approximations, 
we set wq equal to the present-day value of the equation 
of state that is found numerically for a = 1). Despite 
the fact that WEfst contains one more fitted parameter 
than wet, it is clear from Fig. 6 that Wfn agrees better 
with the numerical behavior of w at high redshift and 
shows comparable accuracy at low redshift. Note, how- 
ever, that if we restrict the redshift range to much lower 
z (e.g., < z < 1), and recalculate w\ to give the best fit 
for tWEfst only within this redshift range, then lOEfst gives 
better agreement than u>fi t , although both expressions 
are accurate to within 2% over this redshift interval. 
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Neither of these results is unexpected. The expres- 
sion for wg t in equation (|36[) is derived by perturbing 
the high-redshift evolution of <f>, so we expect it to be 
more accurate at higher redshifts than u>Efst- On the 
other hand, the expression for WEfst was derived from an 
empirical fit at low z, and so it is no surprise that this 
expression can be made to work better at low redshifts. 

While we have tested our results against a fiducial flat 
model with f2 m o = 0.3, nothing in our analysis depends 
on this assumption. For obvious reasons, the agreement 
between our analytic results and the true evolution of 
Ptj, and w will be better than the ones presented here if 
Sl m o > 0.3, and worse if fl m o < 0.3. 

Although we have concentrated specifically on the in- 
verse power-law potentials in this paper, it is obvious 
that the techniques exploited here can be applied to any 
quintessence model in which the quintessence energy den- 
sity is subdominant at early times and dominates at low 
redshift. 



FIG. 6: The ratios WEtst/w (black, solid), and wat/w (red, 
dotted) for the a — 1 case. 
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